Flux qubit as a sensor for a magnetometer with quantum limited sensitivity 
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We propose to use the quantum properties of a superconducting flux qubit in the construction 
of a magnetometer with quantum limited sensitivity. The main advantage of a flux qubit is that 
its noise is rather low, and its transfer functions relative to the measured flux can be made to 
be about 10mV/$o, which is an order of magnitude more than the best value for a conventional 
SQUID magnetometer. We analyze here the voltage-to-flux, the phase-to-flux transfer functions and 
the main noise sources. We show that the experimental characteristics of a flux qubit, obtained in 
recent experiments, allow the use of a flux qubit as magnetometer with energy resolution close to 
the Planck constant. 

PACS numbers: 74.50.+r, 84.37.+q, 03.67.-a 



Josephson-junction qubits are known to be candidates 
for solid-state quantum computing circuits^. However, 
owing to their unique quantum properties these devices 
undoubtedly can be used as sensitive detectors of dif- 
ferent physical quantities, such as quantum environmen- 
tal noised or low frequency fluctuations of the junction 
critical current*. Here we propose to use a Josephson- 
junction flux qubit as a sensitive detector of magnetic 
flujii We show that the present state-of-art allows one 
to obtain the energy sensitivity of such a detector in the 
order of the Planck constant. 

A flux qubil£*§i! consists of three Josephson junctions 
in a loop with very small inductance L, typically in the 
pH range. This ensures an effective decoupling from the 
environment. Two junctions have an equal critical cur- 
rent I c and (effective) capacitance C, while those of the 
third junction are slightly smaller: al c and aC, with 
0.5 < a < 1. At sufficiently low temperatures (typically 
T « (10 ~ 30) mK) when $ x , the external flux applied to 
the qubit loop, is in the close vicinity of $o/2 ($o = h/2e 
is the flux quantum, h is the Planck constant) the system 
has two low- lying quantum states E— and E + . The en- 
ergy gap of the flux qubits, A/h = (E + —E-)/h is of the 
order of several GHz. Below we assume ksT « A (fcs 
is the Boltzmann constant), so that the qubit is definitely 
in its ground state . 

For experimental characterization the flux qubit is in- 
ductively coupled through a mutual inductance M to an 
LC tank circuit with known inductance Lt, capacitance 
Ct, and quality Q (Fig-QJ. 
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The resonant characteristics of the tank circuit (fre- 
quency, phase shift, etc.) are sensitive to the qubit in- 
ductance and therefore to the external flux <£v 

It was shown inS that the amplitude v and the phase 
X of the output signal V(t) — vcos(ujt + x) are coupled 
by the equations: 



1 + 4Q 2 C 2 (/ X )) =lW T L%Q 2 
tan X = 2QC(/x) , 
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where Iq is the amplitude of the driving current 
(Ib(t) = lo cos w£), and to and wt are a driving frequency 
and the resonant frequency of the tank circuit, respec- 
tively. 

It is worth noting that the the scheme in Fig. ^ and 
Eqs. ||TJ, (J2J are quite similar to those for a conventional 
RF SQUID. The only difference is in the expression for a 
flux-dependent frequency detuning This depends 

on the qubit parameters aA 
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£0 = (wt-w)/wt, fx = *x/$o - I, and 
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FIG. 1: Flux qubit coupled to a tank. 



with 77 = 2£jA/A and 7 = MI Q/<&q. The expression 
for A, which depends on a, Ic and C is given inS. 

Therefore, the main effect of the qubit-tank interac- 
tion is a shift of the tank resonance. This results in a 
dip in the voltage-to-flux and phase-to-flux characteris- 
tics which have been confirmed by experiments^. 

Theoretical phase-to-flux x(fx) (PFC) and voltage-to- 
flux (VFC) v{f x ) dependencies at resonance w — ujt, are 
shown in Fig. [21 for three values of the amplitude of the 
bias current Iq. The graphs hasve been calculated from 
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FIG. 2: (Color online). Tank phase X (a) and tank voltage FIG. 3: (Color online). Phase-to-flux x * (a) and voltage-to- 
V (b) vs bias flux f x for three values of bias current I . The fl ux (b) transfer functions for three values of bias current 



gap frequency A/h = 2GHz. 



Iq. The gap frequency A/h = 2GHz. 



Eqs.QJ, (J2J) for the following qubit-tank parameters: I c = 
400 nA, a = 0.8, L = 40 pH, L T = 50 nH, Q = 10000, 
w t /2tt = 100 MHz, A/h = 2GHz, and k = 10~ 2 . 

The advantage of a qubit magnetometer over a conven- 
tional SQUID magnetometer is in the very steep depen- 
dence of its VFC and PFC. In the flux locked loop opera- 
tion of a magnetometer the working point is set at a fixed 
value of <&x where the slope of VFC or PFC is maximum. 
The output signal 8V is proportional to the change 6&x 
of the measured flux. In principle two modes of detection 
are possible: voltage mode, where 5V — VqSQx, and the 
phase mode, where SV = x$<5$x- The qubit transfer 
functions \<s> = vd\/d$x and V$ = dv/d^x are shown 
in Fig. |3| for the same qubit-tank parameters as those 
used in Fig. [21 It is seen that qubit transfer functions 
can exceed 10 mV/$o- This value should be compared 
with 1 mV/$o, the best value obtained for a DC SQUID 
with additional positive feedback^. 

The flux and energy sensitivity depend on the main 
noise sources, which come from the low frequency fluc- 
tuations of the junction critical current, Si c , and from 
the voltage noise, Sy and the current noise, Si of the 
preamplifier, where Si c , Sy and Si are the correspond- 
ing spectral densities. The fluctuations of Ic result in 

1 /2 1/2 

the fluctuating flux in the qubit loop S ^ r = LSj . For 



I c = 400 nA, junction area A = 0.12/im 2 , T = 0.1 K 
we estimate for three-junction flux qubit (see Eq.18 in 3 ) 
Sl /2 Ic » 2 x lO^o/Hz 1 / 2 at 1 Hz. We will see that 
this value is almost one order of magnitude smaller than 
the noise from a preamplifier. Therefore, the self noise 
of the qubit can be neglected. The contribution of the 
voltage noise of the preamplifier to the flux resolution 
referred to the input is S]£ = Sy/V£ or S]£ = Sy/x% 
depending on the detection mode. The current noise of 
preamplifier which is related to its noise temperature Tjv, 
5/ = AMbTn/Rt, {Rt — utLtQ), contributes via two 
mechanisms. The first one comes from magnetic coupling 
between the tank inductance and the inductance of the 
qubit loop S^, = M 2 Q 2 Si. This contribution cannot be 
separated from the measured flux. The second mecha- 
nism contributes through a voltage noise induced by the 
current noise of the preamplifier across the dynamic resis- 
tance of the tank Sg = R 2 D Si/Vg (or Sg = R 2 D Si/x%), 
where Rd = dv/dlo- By combining these three mecha- 
nisms we obtain for the flux sensitivity: 

5$ = M 2 Q 2 Si + S v /Vi + R 2 D Si/Vi (5) 

where Rd is approximately equal to Rt, the resistance of 
unloaded tank. In the case of the phase mode detection 
we should substitute in JSJ) for V$. 
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FIG. 4: (Color online). Flux sensitivity of qubit magnetome- 
ter. Phase detection mode is shown by boxes and circles. 
The stars and triangles are for the voltage detection. Two 
straight lines show the level of S$,i for k = 0.01 (upper line) 
and k — 0.005 (lower line). 

For the estimation we take S\f 2 = 0.2 nV/Hz 1 / 2 , and 
Tjv = O.lKii, with the other parameters being the same 
as for Figs. [3 El The inspection of Eq. JSJ shows that 
the main contribution to the flux noise comes from the 
first term: 

S„ r = M 2 Q 2 Sj = k 2 Q 4kBTNL (6) 
This contribution does not depend on the position of the 

1/2 

working point and for our parameters it gives = 
2.8 x 10~ 7 $o/Hz 1/2 . The influence of the last two terms 



in Eq. (J5J depends on the position of the working point 
and on the bias current amplitude Iq. In general, the 
contribution of these terms is nonnegligible. The total 
flux noise dependence on the amplitude of bias current is 
shown in Fig. 0] As seen from the figure phase detection 
in general is more favorable than voltage detection. It 
gives lower noise and is weakly sensitive to Iq- The flux 
resolution can be improved by a decrease of S$.i (see Eq. 
©) upon optimization of k, Q, ujt, or L. However, it 
does not necessarily lead to a decrease of the total noise, 
since the transfer functions also depend on these param- 
eters. (See, for example, the k — 005 curve in Fig. 0] 
for voltage detection) . An increase of the bias frequency 
lot can also give an improved flux resolution. However, 
for the qubit to remain in the ground state the condition 
lot « A//i should hold. We also made calculations 
for ujt = 200MHz with k=0.01, Q=1000, with other pa- 
rameters being unchanged. We obtain at Jo = 200pA 
for the phase detection S^ 2 = 1.6 x 10" 7 $ /Hz 1 / 2 , 
which for L = 40pH corresponds to the energy sensi- 
tivity s = S&/2L = 1.3 x 10 _33 J/Hz=2/i. These values 
should be compared with those for conventional SQUIDs: 
4 /2 « lO-^o/HzVa, £ « lO-^J/Hai 3 .. 

In summary, we have shown that a superconducting 
flux qubit can be developed as a sensor of magnetic flux 
with an energy sensitivity close to the Planck constant. 
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